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We analyze a controllable generation of maximally entangled mixed states of a circuit 
containing two-coupled superconducting charge qubits. Each qubit is based on a Cooper 
pair box connected to a reservoir electrode through a Josephson junction. Illustrative 
variational calculations were performed to demonstrate the effect on the two-qubits en- 
tanglement. At sufficiently deviation between the Josephson energies of the qubits and/or 
strong coupling regime, maximally entangled mixed states at certain instances of time is 
synthesized. We show that entanglement has an interesting subsequent time evolution, 
including the sudden death effect. This enables us to completely characterize the phe- 
nomenon of entanglement sharing in the coupling of two superconducting charge qubits, 
a system of both theoretical and experimental interest. 



1 Introduction 

There have been remarkable advances in the quest to build a superconductor-based quan- 
tum information processor in recent years [U |2] and one of the greatest scientific and 
engineering challenges of this decade is the realization of a quantum computer. In this 
context, a solid-state system is highly desirable because of its compactness, scalability and 
compatibility with existing semiconductor technology. One of the physical realizations of a 
solid-state qubit is provided by a Cooper pair box which is a small superconducting island 
connected to a large superconducting electrode, a reservoir, through a Josephson junction 
[3]. Superconducting charge qubits (Cooper pair boxes) are a promising technology for 
the realization of quantum computation on a large scale [H El [6l [TJ [8] . 

Using simultaneous measurement and state tomography, entanglement between two 
solid-state qubits has been demonstrated p. The results demonstrate a high degree of 
unitary control of the system, indicating that larger implementations are within reach. 
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These results are promising for future solid-state quantum computing. For conventional 
fault-tolerant quantum computing, the quantum states should have a high level of purity, 
preferably being as close to a pure state as possible. When the qubit is coupled to an 
environment it is subject to decoherence, which will typically result in a completely mixed 
state [lOj. However, a qubit initially in a completely mixed state can be purified by mea- 
surement. Therefore, it is desirable on both fundamental and practical grounds to study 
maximin entangled state generation and entanglement dynamics in a time-dependent 
sense. One of the next major steps towards building a Josephson junction quantum com- 
puter prototype will be the demonstration of controllable coupling between the qubits. 
At this end, it seems that a quantitative link between the degree of disentanglement and 
the amount of the energy transferred between the system of interest and its environment 
is still missing. Investigation of entanglement control in such systems would therefore be 
an important contribution to the present suite of experimental controls. 

In recent years quantum entanglement has found many exciting applications that 
have considerable bearing on the emerging fields of quantum information and quantum 
computing [HI [T21 [131 [H Moreover, besides this fundamental aspect, the interest 
in entangled states has been recently renewed because their properties lie at the heart of 
many potential applications. The generation and reconstruction of quantum states were 
extensively studied in the past theoretically and experimentally [Ml EZl [ISl [IHl [SD]- More 
fundamentally, decoherence processes due to the interaction with internal or external 
noises and entanglement decay in a time- dependent sense have been studied in many 
distinct cases [El [2ll |23l [2l|25l [261 [27] . More recently, Almeida et al. [H] have devised 
an elegantly clean way to confirm the existence of so-called entanglement sudden death, 
that is, entanglement terminates completely after a finite interval, without a smoothly 
diminishing long-time tail. 

This paper examines the generation of maximally entangled mixed state of two-coupled 
Josephson charge qubits using a common pulse gate. We present various examples in 
order to monitor different regimes of synthizing the maximally entangled mixed states and 
entanglement dynamics. In principle, by proper adjustment of the initial state parameters, 
we can always find suitable values of characteristic energies of the Cooper pairs and 
coupling energy which can be used to suppress the decay of entanglement. This analysis 
is carried out with generalized time-dependent density matrix, yielding a generalized 
dynamical two-qubit model. As we are going to show, we may take this advantage to 
increase time intervals for maximally entangled states caused by the strong coupling 
regime. 

This paper is organized as follows: In Sec. 2, we will describe the Hamiltonian of 



the system of interest, and obtain the exphcit analytical solution of the master equation 
describing the dynamics of two qubits in the presence of phase decoherence. In Sec. 3, 
by calculating the occupation probabilities of the two qubits, we show that it is possible 
to generate the maximally entangled mixed states of the system in different situations. 
In Sec. 4, we discuss the entanglement of the system by virtue of the concurrence in the 
absence or presence of the decoherence. Finally, Sec. 5 presents the conclusions and an 
outlook. 
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Figure 1: Illustration of two capacitively coupled Josephson charge qubits. The circuite 
consists of two charge qubits that are coupled by an on-chip capacitor Cm O [28] . 

2 Two coupled charge qubits 

Here, we briefly discuss the general formalism to characterize the dynamics of two-coupled 
superconducting charge qubits (Cooper pair boxes connected to a reservoir electrode 
through a Josephson junction). For a more detailed discussion we refer the reader to 
Ref . P, |30] . We consider two charge qubits and couple them by means of a miniature 
on-chip capacitor. The read-out of each qubit, in this case, is done similar to the single 
qubit read-out and connect a probe electrode to each qubit. External controls that we 
have in the circuit are the dc probe voltages and V^^, dc gate voltages Vg^ and Vg^, 
and pulse gate voltage Vp (see figure 1). The information on the final states of the qubits 
after manipulation comes from the pulse-induced currents measured in the probes. By 
doing routine current- voltage-gate voltage measurements, we can estimate the capaci- 
tances. We then perform state manipulation and demonstrate qubit-qubit interaction. 



The Hamiltonian of the system in the charge representation can be written as 

oo oo p p 

H = fiJ2Yl Viini,n2)Su - -f- {s,2 + Su) - -f- (Srs + ^24) , (1) 

ni=0 n2=0 

where Sn = |ni, nz) (ni, na], 6^12 = \ni,n2){ni + l,n2\, 6^13 = |ni, n2)(ni, ^2 + 1|, 5^34 = 
|?T,i, n2 + 1) {rii + l, n2 + l\, and 5*24 = + 1, ^2) (rii + l, n2 + l|. The parameter rii{ni,n2) = 
Eci {ng^ - rii)^ + {rig^ - ^2)^ + Em {ug^ - rii) {rig^ - ^2) . Here, rii and ^2 (rii, ^2 = 
0, ±1, ±2, ...) are the numbers of excess Cooper pairs in the first and the second Cooper 
pair boxes, and Ug-^^ = {Cg-^^Vg^^ + CpVp)/2e are the normahzed charges induced on the 
corresponding qubit by the dc and pulse gate electrodes. The eigenenergies, Ek {k = 
0,1,2,...), of the Hamiltonian form 2e-periodic energy bands corresponding to the 
ground {k = 0), first excited {k = 1), etc. states of the system. Ec^^, E^.^ and E^. give the 
characteristic energies of Cooper pair of the first qubit. Cooper pair charging energy of 
the second qubit and the coupling energy, respectively. 
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where C^^ ^ are the sum of all capacitances connected to the corresponding Cooper pair 
box including the coupling capacitance Cm and e is the electron charge. 

If the circuit is fabricated to have the following relation between the characteristic en- 
ergies: Ej^ 2 ~ < 2, then one can use a four-level approximation for the description 
of the system (|00), |01), |01) and |11)) around Ug-^ = Ug^ = 0.5 while other charge states 
are separated by large energy gaps. In this basis, the two charge qubits system behaves 
as a single four-level system which can be used as a new basis for the Hamiltonian ([1]). 

The time evolution of the system density operator p{t) can be written as [311 [321 [33] 

= (3) 

where 7 is the phase decoherence rate. Equation ([3|) reduces to the ordinary von Neumann 
equation for the density operator in the limit 7 — 0. The equation with the similar form 
has been proposed to describe the intrinsic decoherence [3l]. Under Markov approxima- 
tions the solution of the master equation can be expressed in terms of Kraus operators 
[35] as follows 

p(t) = exp (^-iH?j exp ("y^') P(0) exp (-y exp (iHt^ H"^ 

= \^ i2!Z_M™(t)p(o)Mt™(t), (4) 
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where p(0) is the density operator of the initial state of the system and M"^ are the Kraus 
operators which completely describe the reduced dynamics of the qubits system, 

= ^™ exp{-iHt) exp (^-y^^^ • (5) 

Equation (jlj) can also be written as 

p{t) = exp(^-iHt^exp(^-^H^^ {e^'''p{0)}exp(^-^H^^exp (6) 
= Pii,zfcWlu)(^^l- 

We define the superoperator S'mp(O) = '-fHp{0)H and choose arbitrary initial state of 
the two charge qubits. The notation = (g) |j)2, is used, where |0)i(2) and |0)i(2) 
are the basis states of the first (second) qubits and Pij^ikif) = {ij\p(t)\lk) corresponds the 
diagonal {ij = Ik) and off-diagonal {ij 7^ Ik) elements of the final state density matrix 
p{t). From here on, for tractability of notation and without loss of generality, we denote 
by Pij(t) = Pij,ijif)i the probability of finding the two-coupled charge qubits in the state 



3 Creation of maximally entangled mixed states 

We now apply the above results to study the time evolution of the occupation probabil- 
ities with different values of the system parameters. In pure-state case, there has been 
considerable debate over the entanglement properties of certain types of states ^llj. In 
this paper, we are interested in the case in which the final state of the coupled charge 
qubits p{t) is a maximally entangled mixed state |19j . 

Nakamura et al. ^36j investigated the temporal behavior of a Cooper-pair box driven 
by a strong microwave field and observed the Rabi oscillations with multi-photon ex- 
changes between the two-level system and the microwave field. Here, the occupation 
probabilities as functions of the scaled time At are schematically shown in Fig. [21 Note 
that the populations of the four states exist but Poo(^) well as Pio(t) oscillate between 
and 1, while Poi(^) ^^'^ Pii(^) oscillate with smaller amplitudes. It should be noted that 
the occupation probabilities results are drastically different when we consider different 
initial state settings. To analyze the effect of the system parameters on the occupation 
probabilities for the present system we consider two different cases. One when the the 
Josephson energies of Cooper pair are different while the second case is the strong inter- 
action regime. This will be seen in figure [3] and figure |H As an example of the creation of 
the two-particle maximally entangled state is shown in Fig. [31 The results of this figure 
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Figure 2: The time evolution of the occupation probabihties Pijit), where Pooit) (sohd 
curve), Poiit) (dashed curve), PiQ(t) (dot-dashed curve) and Pii(t) (dotted curve). The 
initial state of the two-charge qubits is assumed to be p(0) = |00)(00| and the parameters 
used in these figures are Eji = Ej2 = 30 and Em = 6. 

are obtained for parameters Eji = 30/ieV, Ej2 = 5/ieV and E^ = 6/ieV. The way to 
determine experimentally the qubits' Josephson energies Eji and Ej2 has been described 
inRef. [30] 
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Figure 3: The same as figure 1, but Eji = SOpeV, Ej2 = 5/ieV and E^ = QpeV. 

If the system starts from excited state, p(0) = |0,0)(0,0|, we see that the occu- 
pation probabilities of the intermediate states tend to zero at any instant of time and 
the Cooper pair system oscillate only between excited and ground states. If the en- 
vironment is switched off, i.e. 7 tends to zero, and using equation ([6]), at some in- 
stant times r ~ -y^, (ra = 1,2,3,...), we can obtain analytically the values of the diag- 
onal and off-diagonal elements of the density matrix as Poo.ool*^) = Pii,ii('^) = 0-5 and 
Poo,iiit) = Pii^oo(^) = C 7^ 0, otherwise Pij^itif) = 0. Which means that, the final state 
takes the form 





p(t) = - |0, 0) (0, 0| + C(|0, 0) (1, 1| + |1, 1) (0, 0|) + - |1, 1) (1, 1| , (7) 

i.e. the final state ([71) becomes a maximally entangled mixed state [19j. This entangled 
state corresponds to the half-probability peaks in figure [3] (C = 0.13). Depending on 
whether Josephson energy of the first qubit is smaller or larger than the Josephson energy 
of the second qubit, the maximally entangled states are created. Similarly, depending 
on the coupling energy the maximally entangled state is characterized as having short or 
long correlation time. Given enough time, the system will therefore reaches a state where 
both excited and ground states have equal occupation probabilities i.e. the coupled-qubit 
system evolves to the maximally entangled mixed state at the times given by 
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Figure 4: The same as figure 1, but Eji = 30/ieV, Ej2 = l/xeV and Em = GOfieV. 

We turn our attention to consider the case in which both characteristic energies have 
different values taking into account the effect of the coupling energy. For this reason 
we have plotted the function Pij{t) against the scaled time Xt in figures (jl]). To make 
a comparison between this case and the previous one we have to take the values of the 
other parameters similar to that of the previous case. In this case and providing the 
characteristic energies Eji = 30/ieV, and Ej2 = GfieV, we find that the function reduces 
its value to be around ~ 1 and 0.5 for the excited state while the ground state prob- 
ability oscillates between and 0.5. Furthermore, we realize there is a long period of 
the interaction time in which the probability of the excited state equals the ground state 
probability, with perfect symmetric fluctuation pattern around 0.5, see figure (jl]). Which 
means that, with these setting, we obtained long lived maximally entangled mixed state 



(in this case ( = 0.19). In the meantime if we exchange the values of the characteristic 
energies, we observe there is no big change in the figure shape, except some decreases in 
the fluctuations number. It is to be noted that the maximally entangled state in this case 
is lived longer compared with significantly short time in the previous case (see figures [3] 
audi]). 

One might now raise the following notes: taking a strong coupling regime where the 
coupling energy is strong enough, one can obtains maximally entangled mixed states at 
some instant times in a periodical manner. But when the coupling energy between the two 
qubits is weak, the period becomes shorter. Also, if the deviation between the Josephson 
energies is substantially large, the maximally entangled states can be generated. It is 
worth noting that when the coupling energy tends to zero, our model becomes similar 
to that of a beam splitter model. In such a case, we cannot obtain maximally entangled 
states. 

The above discussion clearly shows that the maximally entangled state generation of 
the two charge qubits depends on both the time evolution, Josephson energies of both 
charge qubits and coupling energy. The considerations of experimental observability of 
the entangling power discussed in [3] are valid in the context of the present work. 

4 Entanglement 

Having established the existence of maximally entangled states p{t), in section 3, now 
we try to answer the following question: how does the entanglement of the two charge 
qubits system evolve? To answer this question, one first needs a formal definition of 
entanglement. Currently a variety of measures are known for quantifying the degree of 
entanglement in a bipartite system [371 EH EHl HOj HH HJl |13]. A convenient measure of 
entanglement for a two-qubit state p(t) is the concurrence Cp (t) , given by 

Cp (t) = max {0, Ai - A2 - A3 - A4} , (9) 

where Ai > A2 > A3 > A4. We denote by Aj the square roots of the eigenvalues of 
p {cTy ® (Ty) (p)* (cTy ® ay) , here ay is the second Pauli matrix and the conjugation occurs 
in the computational basis (|00), |01), |10), |11)). Cp (t) quantifies the amount of quantum 
correlation that is present in the system and can assume values between (only classical 
correlations) and 1 (maximal entanglement). 

In figure [5l we plot the concurrence as a function of the scaled time assuming that the 
two-coupled superconducting charge qubits start from their excited states. The maximum 
value of the entanglement decreases as the coupling energy is decreased. As time goes on. 




Figure 5: The time evolution of the concurrence as a function of the scaled time Xt. The 
parameters are E jj=30fieV , Ej2=2fieV, 7 = and different values of the coupling 
energy, where Em=200fieV (dotted curve), Em=60 ^eV (dashed curve) and Em=5^eV 
(solid curve). Regions of the entanglement sudden death are painted in gray. 



the entanglement reaches zero value in a periodic way, this period decreases as the coupling 
energy decreases. In the uncoupled situation each qubit oscillates with its own frequency 
and Cp (t) =0. It is interesting to note that, the maximum entanglement is achieved 
at specific choices of the interaction time i.e. the entanglement content corresponding 
to specific choices of the interaction time and large values of the coupling energy. We 
therefore consider the question of how the coupling energy affect the entanglement of 
system. In relation to that discussion, it is useful to examine the effect of the characteristic 
energies with fixing the coupling energy of the two charge qubits that change their state 
during the transition. That criterion is related to, but clearly distinct from, the question 
of quantifying how maximum a quantum state is. Also, as can be seen from the graphs 
figures [3] and 0, the entanglement vanishes at the time at which the population of the 
symmetric state is maximal. 
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Figure 6: The same as figure 4 but E j2=5 fieV. 

Put differently, with a small difference between the characteristic energies of the 
Cooper pairs (Josephson energy of the second qubit, E j2=5^eV) and still kept the pa- 
rameter values of figure [5l we show that the entanglement features are visibly worsened, 
as in this regime oscillations are faster than the previous case (see figure [6]) in which big 
difference between the characteristic energies of both qubits are considered. One has still 
zero entanglement due to the time development, but with a short period of the interaction 



time which is roughly given by Xt ~ (n = 1,2,3, ...). Similar to the previous case, 
one will have a very small amount of entanglement when the coupling energy decreases 
and this amount disappear completely when the coupling energy decreases further. In- 
deed, the comparison of plots figure O and figure El demonstrates that the entanglement 
in both cases has somewhat similar behavior corresponding to different values of coupling 
energy. The deviation value of the Josephson energies of the Cooper pairs effect on the 
entanglement is particularly pronounced as this deviation is much bigger. 

Conventionally, maximally entangled state emerge from the coupling of the two qubits 
by a small island overlapping both Cooper pair boxes, i.e. two-coupled superconducting 
charge qubits. As such, the phenomenon is the result of many-particle dynamics, often 
described by a simple interaction model. In the single-qubit case novel features appear 
which are due to the coherent microscopic dynamics. Our study allows us to identify 
the dependence of these features on different parameters of the system, thereby giving us 
insight into how maximally entangled state and hence maximum entanglement arise from 
dynamics in this particular coupling process. 

The quantum features of many systems decay uniformly as the result of decoherence 
and much effort has been directed to extend the coherence time of these qubits. However, 
it has been shown that under particular circumstances where there is even only a partial 
loss of coherence of each qubit, entanglement can be suddenly and completely lost [131 HI] . 
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Figure 7: The same as figure 5 but Em=200fi and different values of the decoherence 
parameter, where 7 = 0.01 (solid curve), 7 = 0.1 (dotted curve) and 7 = 0.8 (dotted- 
dashed curve). 

This has motivated us to consider the question of how decoherence effects the scale 
of entanglement in the present system. The decoherence time due to the coupling to the 
vacuum via the probe junction can be estimated to be roughly 100 ns at the resonant 
condition. For the most of the experimental devices [7j, this probing time restricted the 



upper limit of the decoherence time, since the probe junction attached to the box had a 
samphng time of typically 8 ns. Once, the environment has been switched on, i.e., 7 7^ 0, 
it is very clear that the decoherence plays a usual role in destroying the entanglement. In 
this case and for different values of the decoherence parameter 7, we can see from figure 
d?]) that after the onset of the interaction the entanglement function increases to reach its 
maximum showing strong entanglement. However its value decreases after a short period 
of the interaction time to reach its minimum. The function starts to increase its value 
again however with lower local maximum values showing a strong decay as time goes on. 
It is interesting to remark that decoherence due to normal decay is often said to be the 
most efficient effect in physics. Which means that, the entanglement increases rapidly, 
then approaches to a minimum value in a periodic manner. Also, from numerical results 
we note that with the increase of the parameter 7, a rapid decrease of the entanglement 
(entanglement sudden death) is shown [35]. 

Since the discovery of entangled sudden death [351 125] , a large number of instances 
of this surprising effect have been identified in the theoretical literature [HI H5l HHl HZ] ■ 
In general, decay takes infinitely long, so one can wait any length of time. However, if 
the entanglement reaches zero in a finite time [35l [25] the game is over. Thereafter no 
distillation process exists that will recover any useful feature of entangled quantum joint 
coherence for use in quantum computing or communication [3B]- These results should 
mark an important consideration in the design and operation of future quantum infor- 
mation networks. Also, properties of such entanglement decay depends on the coupling 
energy and qubits Josephson energies. Even if achieving the maximally entangled state is 
not possible in the presence of the decoherence, one can argue that finding a maximally 
entangled state can, under certain conditions, be achieved for a short interaction time. It 
is also worth noting here that we have used the simple model of a two-Cooper pair box 
problem, which represents a proper physical system that can be used as a qubit. This 
system is of a great interest because it offers the possibility of scaling to a large number 
of interacting qubits. 

5 Conclusion 

In the important context of quantum state engineering and characterization, we have stud- 
ied the entanglement properties of the special class of two-coupled superconducting charge 
qubits. Clear physical interpretations for the maximally entangled state generation and 
entanglement found for certain parameter regimes of the system have been provided. In 
a strong coupling regime, two different types of maximally entangled states of the system 



have been created. This helps the comprehension of the quantitative results achieved by 
the use of the concurrence. Moreover, we find peculiar entanglement characteristics which 
are unique to this system, and which we trace back to the interplay of the various time 
scales of the dynamics. Our results suggest that the maximally entangled mixed state is 
generic in different situations for the coupled-charge qubits system, and that developing 
entanglement theory under other sorts of restrictions is a promising direction for further 
study. In a more general context, our results provide further insight into the coupled 
dynamics of superconducting charge qubits in the spirit of the experimental realization 
[31 [TJ [29] and may provide a useful maximally entangled states source in the exploration 
of various quantum-information processing. 
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